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Abstract 

The phenomenology of velocity statistics in turbulent flows, up to now, relates to 
different models dealing with either signed or unsigned longitudinal velocity incre- 
ments, with either inertial or dissipative fluctuations. In this paper, we are concerned 
with the complete probability density function (PDF) of signed longitudinal incre- 
ments at all scales. First, we focus on the symmetric part of the PDFs, taking into 
account the observed departure from scale invariance induced by dissipation effects. 
The analysis is then extended to the asymmetric part of the PDFs, with the spe- 
cific goal to predict the skewness of the velocity derivatives. It opens the route to 
the complete description of all measurable quantities, for any Reynolds number, 
and various experimental conditions. This description is based on a single universal 
parameter function T>(h) and a universal constant 1Z*. 
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1 Introduction 

In the field of turbulence, a significant effort has been devoted to the anal- 
ysis of the scaling behavior of structure functions {(5iu) q ), where 5gu is the 
longitudinal velocity increment between two points separated by a variable 
distance £ [1] . However, a better strategy may be to concentrate on the prob- 
ability density functions (PDFs) of 5eu, rather than on a set of moments [2]. 
Accordingly, this work deals with the description of the PDFs of 5eii, where 
the scale £ spans the entire range of excited scales of motion (from the integral 
far down to the dissipative scales). Experimental and numerical observations 
have provided the evidence that the PDFs of 5eu are increasingly stretched 
as £ decreases, while they are almost Gaussian at the large scales where the 
turbulence is stired [1]. This feature is known as intermittency. Moreover, 
Chevillard et al. [3] have recently argued that this stretching is largely en- 
hanced in the near-dissipation range, leading to extremely high fluctuation 
level for the velocity gradients. Another essential feature lies in the significant 
asymmetry, or skewness, of the PDFs. This skewness to be non-zero is heuris- 
tically connected to the vortex folding and stretching (irreversible) process, 
which drains energy from the large to the small scales, and hence, plays a 
central role in turbulence. From a theoretical viewpoint, a quantitative de- 
scription of the skewness is still missing. In this context, our motivation is 
to present a synthetic description of the PDFs of Sgu, which encompasses the 
combined effects of intermittency and skewness. To do so, the PDF of 5eu is 
split into a symmetric V^ u and an asymmetric Vj lU part. First, the focus is on 
Vt u , which also represents the PDF of the magnitude of 5(U. We show that 
Vf lU is suitably described by a multifractal picture of turbulence dynamics [1], 
which incorporates finite-Reynolds-number effects. The analysis is then ex- 
tended to Vj lU with the specific goal to describe the skewness phenomenon via 
a quantitative estimate of the skewness factor (5eu 3 ) / {S(U 2 ) 3 ^ 2 as a function 
oil 



2 Statistics of longitudinal velocity increment magnitude: Model- 
ing the symmetric part of the Probability Density Function 

From a general point of view, the PDF of the longitudinal velocity increments 
V$ eU can be split into a symmetric V^ lU (i.e. even) function and an asymmetric 
Vj tU (i.e. odd) function in the following way: 

= Vt tu {5,u) + Vi eu (5 e u) . (1) 

The PDF of the longitudinal velocity increment magnitude V\s eU \(\5tu\) = 
'Ps t u(\de u \) +T ? s e u(~\^e u \) — 27-^(1(5^1) shows that the symmetric part of the 
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PDF of 5eu describes the magnitude statistics. Notice that neither Vf lU nor 
VJ U can be interpreted as a PDF of a random variable. 

Let us first focus on the symmetric part of the PDFs of the longitudinal ve- 
locity increments Vf{5^u). In the inertial range, the multifractal formalism 
[4], which a priori pertains in the limit of infinite Reynolds number, states 
that velocity is everywhere singular, the longitudinal velocity increments at 
scale I behaving locally as £ h , where the Holder exponent h takes value in a 
finite interval [/i mm , Vax]- When the Reynolds number lZ e = oLjv is finite 
[L is the correlation length scale, a = \J ((5lu) 2 ) and v the kinematic vis- 
cosity), Paladin and Vulpiani [5] have argued that the dissipative scale, that 
is supposed to separate the inertial and the dissipation scaling ranges, is not 
unique in the presence of intermittency and is likely to depend on h. Using 
these arguments, Nelkin [6] predicted the moments of velocity gradients, i.e. 
((d x u) q ). The phenomenological consequences on the energy power spectrum 
were studied by Frisch and Vergassola [7] who proved the existence of an inter- 
mediate dissipative range. Meneveau [8] further investigated the behavior of 
the structure functions in that transitory range. Recently, Chevillard et al. [3] 
revisited the behavior of longitudinal velocity increments in the intermediate 
dissipative range and showed, among other new predictions, that the width 
[?7_, of this range of scales behaves non trivially with the Reynolds number, 
i.e. ln(77 + /V) ~ y/\nTZ e . 

To provide a complete statistical description of longitudinal velocity incre- 
ments statistics, one needs to model the probability law of the stochastic vari- 
able \8(u\. As originally proposed by Castaing et al. [9], within the propagator 
approach, velocity increments magnitude can be considered as the product of 
two independent random variables, \8(.u\ — x \S\ (in law), where 5 is a zero- 
mean gaussian random variable of variance a 2 and a positive random vari- 
able (see [3] for details). In the inertial range, i.e. rj(h) C ^ C I, where r)(h) is 
the fluctuating dissipative scale, (3e(h) = {£/ L) h can be expressed as a function 
of the singularity strength h that fluctuates from point to point according to 
the probability law Ve(h) ~ (£ / 'L) 1 ~ T>( - h \ Note that the exponent h and the 
parameter function V{h) gain the mathematical status of Holder exponent 
and singularity spectrum in the inviscid limit (72 e — > +oo). The dissipative 
scale r}{h) fluctuates according to : r](h) = L(7£ e /7?*)~ 1 / ( ^ +1 - ) , where the con- 
stant 72* is necessary to be consistent with experimental and numerical data 
[10,11,12]. More precisely, in a monofractal description of velocity fluctuations 
(Kolmogorov K41 theory [1]), one can show [3] that the Kolmogorov constant 
c K = ((5 e u) 2 )/(e) 2 / 3 f/ 3 = (ft*/15) 2/3 . Actually in this simplified monofractal 
framework, the average local dissipation (e) = 15u((d x u) 2 ) = (15/72.*) x a 3 /L 
can be related to the ratio lZ e /(lZ\) 2 = 4/72* (where 1Z\ is the Taylor based 
Reynolds number). We will see in the following that the data are compatible 
with the universal value 72* = 52 (in the presence of intermittency). 
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For scales £ < i](h), the velocity is smooth and Taylor's development applies, 
i.e. 8(u(x) = £d x u(x). In the multifractal description, using a simple continuity 
argument with the inertial range behavior [6] yields f3g(h) = (£ / 'L)(r)(h) j 'L) 11 ^ 1 
and Ve(h) ~ (r](h) / L) 1 ^^. Then, we impose that the function (3 e be con- 
tinuous and differentiable at the transition, following a strategy already used 
in a slightly different form in Ref. [8], and which is inspired from an elegant 
interpolation formula originally proposed by Batchelor [13], independently de- 
rived in a field theoretic approach [14]. In this framework, a single function 
(3e(h,1Z e /Tl*) covers the entire range of scale: 



(3 e (h,n e /n* 



ay 



21 (1-70/2 



(2) 



and 



e \l-V(h) 



V e {h,K e /n*,V) = 



tt) 



Z{t) 



(V(h)-l)/2 



(3) 



h(iy)" 

where Z(£) is a normalization factor such that J*^™** Vg(h, 1Z e /1Z*, V)dh = I. 

From Eqs. (2) and (3), one can derive analytical predictions for the moments 
of the longitudinal velocity increment modulus, i.e. (I^wl 9 ), the energy power 
spectrum (which is linked to the Fourier transform of the second order moment 
[15,16]) and the symmetric part of the longitudinal velocity increments PDF. 
This approach has also been successfully applied in the Lagrangian framework 
in which the PDFs are symmetric [17,18]. 

As advocated in Ref. [19], the magnitude cumulant analysis provides a more 
reliable alternative to the structure function method. The relationship between 
the moments of \8gu\ and the cumulants C n (£) of In \8eu\ reads 



(M p )=exp E^nW" 



(4) 



\n=l 



In Fig. 1, we report the results of the computation of the magnitude cumulants 
C n (£) of various experimental velocity signals. Actually we have plotted the cu- 
mulants of In instead of In \ 8iu\, so that they vanish at the correlation length 
scale L (as the signature of Gaussian statistics). When both the cumulants and 
\n(£/L) are renormalized by \n(lZ e /lZ*) : all the curves collapse on a universal 
linear function in the inertial range (when \n(£/ L) / \n(1Z e /1Z*) > —3/4, see [3]) 
of slope c n . Let us notice that in this representation, for any Reynolds num- 
ber, \n(rj K /L)/\n(TZ e /Tl*) = -3/4 and ln(A/L)/ \n(R, e /K*) = -1/2, where 
t)k and A are respectively the Kolmogorov and the Taylor scales. For the 
first order cumulant (Fig. 1(a)), disregarding large scale anisotropy leading 
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Fig. 1. Magnitude cumulant analysis of various experimental longitudi- 
nal velocity profiles : (o) Turbulent low temperature gaseous helium jet 
for Tlx = 208 [20]; (o) Air Jet for TZ X = 380 [21]; (v) Wind tun- 
nel for TZ X = 2500 [22]. (a) d{£) = (In fa) = (ln\5 e u\) - (ln|<5|). (b) 
C 2 (£) = Var(ln/3 £ ) = (In 2 (3 e ) - (ln/3 £ ) 2 = Var(ln|^u|) " Var(ln|<5|). (c) 
C3 = (In 3 fa) — 3(ln 2 Pi) (In Pi) + 2 (In Pe) 3 . The solid curves correspond to our theo- 
retical predictions (see text). 



to nonuniversal effects, C\ is found very close to 1/3, consistently with K41 
theory [1]. For the second-order cumulant (Fig. 1(b)), the intermittency coef- 
ficient c 2 = 0.025 ± 0.003 is found universal, i.e. independent of the Reynolds 
number and of the experimental configuration. For the third one, C3 can- 
not be claimed to be different from zero (especially at high Reynolds num- 
ber) confirming that the statistics of longitudinal velocity increments are 
likely to be log-normal [19,23]. In the intermediate dissipative range (i.e. 
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-1.1 < la(£/L)/ln(K e /K*) < -3/4), d(£) crosses over towards trivial scal- 
ing; the straight line of slope unity observed at smaller scales means that 
velocity increments become proportional to the scale (Taylor development). 
The behavior of the second-order cumulant is much more interesting and has 
been widely studied in Ref. [3] : a non trivial Reynolds dependent rapid in- 
crease occurs in the intermediate dissipative range, the larger the Reynolds 
number, the more "rapid" the increase. Finally, C^)/ In (lZ e /1Z*) tends to- 
ward a universal value in the far-dissipative range. Note that C 3 (£) displays 
similar behavior. In Fig. 1 are also represented our theoretical predictions 
obtained from the computation of the moments of Infy using Eqs. (2) and 
(3), i.e. ((\n(3 e ) n ) = f^(\n [3 e ) n V e (h)dh. We have used the following set of 
parameters: TV = 52 and a universal log-normal parabolic T>{h) function, 



with c 2 = 0.025 andci = l/3+3c 2 /2 « 0.37 to ensure that Cs = 3^-9^/2 = 1 
in the inviscid limit [1]. The integration limits h niin and h nmx are respectively 
the minimal and maximal values such that T>(h) > 0. Using Eq. (5), we get 
ft-min — Ci — \^2~C2 ~ 0.15 and h max — C\ + V2c 2 ~ 0.59. The different curves 
so-obtained superimpose remarkably well with the corresponding data for the 
first two cumulants which demonstrates that our multifractal description ac- 
counts quantitatively well for the departure from scaling in the intermediate 
dissipative range. Finite Reynolds number effects [23], statistical convergence 
and lognormal approximation can explain some discrepancies between our 
theoretical prediction and the behavior of the third-order cumulant. 



3 Multifractal prediction of the skewness of longitudinal velocity 
increments 

Let us now investigate the statistics of signed longitudinal velocity increments 
through the two estimators: (i) the skewness S(£) = — ((5eu) 3 )/((5iu) 2 ) 3 ^ 2 
and (ii) the asymmetry factor A(£) = — ((5eu) 3 ) / {\Sgu\ 3 ) . The experimental 
estimates of S(£) and A(£) are shown in Fig. 2 in a semi-logarithmic repre- 
sentation. Interestingly, A(£) displays a plateau at about 0.14 in the inertial 
range, whereas the skewness behaves approximatively as a power law. In the 
intermediate dissipative range, both estimators undergo a rapid acceleration, 
very much like what was observed for C 2 (^) m Fig- 1(b)- From a theoretical 
point of view, the third-order structure function is solution of the Karman- 
Howarth-Kolmogorov equation [24] 



V(h) = 1 



(h- Cl ) 2 
2c 2 



(5) 




d((5 e u) 2 ) 
d£ 



(6) 



6 




ln(l/L)/ln(Re/R*) 

Fig. 2. Asymmetry factor A(£) and Skewness S(£) estimated from various experi- 
mental velocity profiles. The symbols have the same meaning as in Fig. 1. The solid 
curves correspond to our theoretical predictions (see text). 

This equation allows us to compute the third-order structure function when 
the second-order one and the average local dissipation are known. A similar 
approach has been performed by Qian in Ref. [25] without taking into account 
the intermittency corrections. 

We have superimposed in Fig. 2 our theoretical predictions to the experimental 
data for S(i) and A(£). Indeed, from Eqs. (2) and (3), one can compute any 
moment of the magnitude of velocity increments (|<^-u| n ) and velocity gradient 
(\d x u\ n ). In particular, we get 

o a 3 15 TZ* 1 f u fn(h)\ 2ih ~ 1)+1 ~ v{h) 
(e) = 15u((d x u) 2 ) = x£--L- IW.) dh , (7) 

W U ' ' LTV ft e Z(0) A min \ L J y 1 

where Z(0) is the limit when I — > of the normalization factor Z{€) appearing 
in Eq. (3). Then, by inserting our description of the second order structure 
function (Eq. (4) for p — 2) together with the prediction of the average local 
dissipation (Eq. (7)) in Eq. (6), we get the third order moment of velocity 
increments ((5^-u) 3 ) at any scale and Reynolds number. As shown in Fig. 2, 
the agreement is very good for distances in between the Kolmogorov and 
Taylor scales r\ K < t < A (without any arbitrary shifts) when using the 
same parameter function V{h) and TV as in our former magnitude cumulant 
analysis. 

Furthermore, we can derive that, when neglecting intermittency corrections 
and setting the viscosity to zero in Eq. (6), we get A(£) = "iV^/TZ* rs 0.145 
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and S(£) = 12/7?* ~ 0.23, in perfect agreement with experimental findings. 
Some discrepancies occur for £ > A, especially for Modane's longitudinal ve- 
locity profile, because (i) of the lack of statistics and (ii) at these scales, one 
has to take into account in Eq. (6) fluctuations of the injection rate of energy 
[26,27]. In the intermediate and far dissipative range, our formalism predicts 
a universal plateau and a Reynolds number dependence for respectively the 
asymmetry factor and the skewness of derivatives, in consistency with Nelkin's 
predictions [6] (see Tab. A.l for precise values). We derive in the appendix A 
the multifractal prediction for the third order moment of the velocity gradient 
{(d x u) 3 ). Experimentally speaking, measuring gradients is still controversial 
mainly because hot wire probe sizes are in general of the order of the Kol- 
mogorov scale [28,29,30,31]. We hope that further experimental studies will 
provide definite test of the validity of our predictions. These preliminary tests 
are nevertheless very satisfactory. 



4 Modeling the asymmetric part of the PDFS 



Let us finally elaborate a formalism to describe the PDF of the signed lon- 
gitudinal velocity increments. To do so, we suggest to model the (signed) 
longitudinal velocity increments in the following way: bgu = ft x A ( (in law), 
where the positive random variable ft is unchanged but A^ is now an inde- 
pendent zero mean random variable of variance a 2 , whose probability a 
priori depends on the scale £. It follows that 

According to the Edgeworth's development [32] , any PDF can be decomposed 
over a "basis" of the successive derivatives of a Gaussian function: 

Pa , W =|A„ M ^(^e-^). (9) 

The symmetric part (even terms) of the PDF of A^, i.e. V£ , is well described 
by a Gaussian 5 noise (as previously stated), which means that Ao(^) = 1 and 
^2n(£) — for n > 1 and every scale £. Furthermore, it can be demonstrated 
that whatever a 2 is, {(Siu) 3 ) = — 6A 3 (£)((ft) 3 ). Hence, A 3 (£) is fully determined 
by the Karman-Howarth-Kolmogorov equation (Eq. (6)). As Eq. (6) is the only 
available constraint on A n , it is quite natural (as a first approximation) to re- 
strict the expansion to A3: A2 n +i(^) = for n > 2 and every scale £. Additional 
statistical equations involving odd moments of 8gu would be needed to give the 
next A 2n +i(^)- This would require further modeling (primarily to get ride of 
pressure terms), which is out of the scope of the present work. Unfortunately, 
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Fig. 3. PDFs of signed longitudinal velocity increments of Air Jet [21] (a-c) and 
Modane [22] (d-f) velocity signals for various scales, (a) and (d) : lnVe(5(u). (b) and 
(e) : (5eu) 3 Ve(5£u). (c) and (f) : (5£u) 4 Ve(5eu). Represented scales (from top to bot- 
tom) : ln(£/L) = -6.0069,-5.3137,-4.6206,-3.9274,-3.2343,-2.5411,-1.8480,0.9246 for 
the Air Jet data and ln(£/L) = -6.4137, -5.6028, -4.6645, -3.6411, -2.7501, -1.8598, 
-0.8685, 0.1226 for Modane data. All curves are arbitrarily vertically shifted for the 
sake of clarity. The solid curves correspond to our theoretical predictions (see text) . 

the previous crude approximation for the odd terms of (9) leads to severe 
pathologies, such as negative probability for rare large events and is not 
consistent with higher-order statistics as Hyperskewness ((<5/u) 5 ) / ((<5/u) 2 ) 5 / 2 
(data not shown). Nevertheless, since the third-order structure function does 
not depend on the precise variance a 2 entering in the third order derivative of 
the Gaussian PDF, we propose to renormalize the variance a 2 of the retained 
odd term (n = 3): a 2 = 0.9a 2 . We thus obtain V Ae (A e ) = V£ e (A e ) + V Ae (A e ), 
where the asymmetric part of the PDF of A e is modeled as 

™ - (t^t • <io) 
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The PDFs of longitudinal velocity increments so-obtained from Eq. (8) are 
shown in Fig. 3 for several scales spanning the inertial and intermediate dis- 
sipative ranges, and compared to the experimental ones for both Air jet Fig. 
3(a-c) and Modane 3(d-f). We see a continuous deformation across scales, 
from Gaussian at the correlation length scale (L) to exponential-like distri- 
butions in the inertial range and ultimately to stretched-exponential when 
dissipation starts acting. This PDF shape evolution is the signature of inter- 
mittency and is remarkably reproduced by our formalism when using the same 
quadratic parameter V(h) function and constant TV = 52 as in Figs. 1 and 2. 
This agreement is emphasized in Figs. (3)(b,d) and (3)(c,f) where respectively 
(Seu) 3 Ve(S£u) and (5eu) 4 Ve(5iu) are represented as a quantitative test of the 
relevance of our formalism to account for the dissymmetrical PDF tails. 



5 Conclusion 

To conclude, we have shown that the evolution across scales of the signed 
longitudinal velocity increments statistics, from the inertial far down the dis- 
sipative ranges, depends only on a universal parameter function T>(h) and a 
universal constant TV that must be seen as a multifractal version of the Kol- 
mogorov constant. In particular, neglecting the intermittency corrections, we 
provide an enhanced phenomenology of turbulence in deriving the value of the 
Skewness S(£) = 12 /TV ~ 0.23 in the inertial range. We have further shown 
that choosing a quadratic form for T>(h) (i.e. the hallmark of an underlying 
lognormal cascading process) provides a very good quantitative description 
of the longitudinal velocity increments PDFs measured in several flows, in 
different geometries and for different Reynolds numbers. This study proposes 
a new formalism, relying on the Edgeworth's development, which opens the 
route to the modeling of velocity increment PDF. New experimental inves- 
tigations of velocity gradients statistics would be welcomed as an additional 
and complementary test of our theoretical multifractal approach. 

We wish to acknowledge P. Flandrin for fruitful discussions. 



A Multifractal prediction of the Skewness of derivatives 

In an infinite domain, or in a finite domain with periodic boundary conditions, 
a Taylor's development of the second order structure function leads to 

((M 2 >H(£d>^ + o(0) )=^muf)--{{dluf)+o{t) . (A.l) 
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Tlx 5(0) .4(0) 



208 


0.35 


0.17 


380 


0.38 


0.17 


2500 


0.50 


0.17 



Table A.l 

Theoretical predictions of the Skewnesses and Asymmetric factors of velocity deriva- 
tives (Eq. (A. 3)) for the three different experiment longitudinal velocity data sets 
presented in Figs. 1 and 2. We have used 1Z* = 52 and a parabolic function for T>(h) 
(Eq. (5)) with parameters C2 = 0.025 and c\ = 1/3 + 3C2/2. 

Inserting the development pointed by Eq. (A.l) into the Karman-Howarth- 
Kolmogorov equation (Eq. (6)) leads to 

((d x uf) = -2u((d 2 x u) 2 ) . (A.2) 

This classical result can be found in Ref . [24] . One may wonder whether our de- 
scription of the second order structure function (Eqs. (2) and (3)) is consistent 
with this development (Eq. (A.l)). In particular, the pre-supposed continuous 
and differentiable transition between the inertial and the dissipative range of 
scale inspired from the work of Batchelor [13] should give a leading term pro- 
portional to £ 3 once inserted in Eq. (6). This property constrains seriously the 
possible form of the transition. The transition form used here benefits of such 
property. We get, with the help of a symbolic calculation software, 

Qua 2 \ 2 r^ax W-2)+i-v(h) 

(A.3) 

where T is a negligible additive term, coming from the Taylor's development 
of the normalization factor Z(£), and given by 

(A.4) 



Eq. (A.3) can be seen as the multifractal prediction of the third order moment 
of the velocity derivatives, using the same transition interpolation form as in 
Eqs. (2) and (3). This is also a prediction for the second order moment of 
the second order derivative of velocity via Eq. (A.2). We gather in Table A.l 
the theoretical values for the Skewness 5(0) = ((d x u) 3 } / ((d x u) 2 ) 3 ^ 2 and the 
Asymmetric factor .4.(0) = ((d x u) 3 ) / ((d x u) 2 ) 3 ^ 2 , i.e. the limit when £ — * of 
our theoretical predictions for the Skewness and Asymmetric factor of velocity 
increments presented in Fig. 2. A specifically designed experiment aimed at 
measuring the fluctuations of longitudinal velocity increments for scale much 



11 



smaller than the Kolmogorov length scale will provide a decisive test of the 
validity of these theoretical predictions. 
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